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Abstract 

In vertebrates, insufficient availability of calcium and inorganic phosphate ions in extracellular 
fluids leads to loss of bone density and neuronal hyper-excitability. To counteract this problem, 
calcium ions are usually present at high concentrations throughout bodily fluids - at concentrations 
exceeding the saturation point. This condition leads to the opposite situation where unwanted 
mineral sedimentation may occur. Remarkably, ectopic or out-of-place sedimentation into soft 
tissues is rare, in spite of the thermodynamic driving factors. This fortunate fact is due to the 
presence of auto-regulatory proteins that are found in abundance in bodily fluids. Yet, many 
important inflammatory disorders such as atherosclerosis and osteoarthritis are associated with 
this undesired calcihcation. Hence, it is important to gain an understanding of the regulatory 
process and the conditions under which it can go awry. In this manuscript, we extend mean-held 
continuum classical nucleation theory of the growth of clusters to encompass surface shielding. 
We use this formulation to study the regulation of sedimentation of calcium phosphate salts in 
biological tissues through the mechanism of post-nuclear shielding of nascent mineral particles by 
binding proteins. We develop a mathematical description of this phenomenon using a countable 
system of hyperbolic partial differential equations. A critical concentration of regulatory protein 
is identihed as a function of the physical parameters that describe the system. 
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I. INTRODUCTION 


In biology, ionic calcium (Ca^^ ) plays many diverse roles including acting as a secondary 
messenger in biochemical cascades and modulating neuronal excitability [1]. Structurally, 
Ca^’*' and inorganic phosphate ions (PO4) are also the main constituents of bones in verte¬ 
brates. For this reason, it is important for organisms to obtain adequate amounts of calcium 
from the environment. 

In normal circumstances, Ca^’*' is plentiful throughout extracellular spaces and in the 
circulatory system, stably existing at concentrations exceeding the saturation point, whereby 
sedimentation is favored | 2 ]. Even under the dangerous condition of hypocalcemia, Ca^"'' 
may still be supersaturated relative to the most thermodynamically stable phase of calcium 
phosphate, hydroxyapatite (HAP), which is the building block of teeth and bones. 

Yet, while the deposition of calcium into bones is desirable, ectopic calcihcation into soft 
tissues is pathological and either causes or exacerbates a variety of inflammatory disorders 
including arteriosclerosis, heart disease, and arthritis [3HS]. So, the regulation of ectopic 
calcium sedimentation is important in maintaining the health of soft tissues. 

When calcium-phosphate solutions are supersaturated, HAP is formed in a multi-step 
process traversing through several intermediate crystalline or psuedo-crystalline states. The 
hrst step in this process is thought to be formation of pre-nucleation clusters [BHH] , which are 
small calcium-phosphate complexes [9]. Although some debate exists [TlElIin], these small 
complexes are thought to be Posner’s clusters (PC), with composition Cag(P 04 )g [H]. A 
combination of experimental and theoretical analyses have conhrmed the stability of PCs [T^ 
[13], their presence in physiological solutions and their consistency with the unit-cell 
structure of calcium-phosphate precipitates P ITl]. For the purposes of this manuscript, we 
will assume that PCs are the fundamental building blocks of larger-scale calcium-phosphate 
clusters and refer to them as monomers. 

These monomers aggregate whereby they nucleate into amorphous spherical post- 
nucleation clusters composed of amorphous calcium phosphate (ACP) [T5IITS] . having a 
calcium to phosphate ratio of approximately three to two [I9]. As long as supersaturation 
persists, and in the absence of regulation, ACP clusters continue to grow by absorbing ad¬ 
ditional monomers into their structure. When ACP clusters become sufficiently large, they 
sediment into the tissue while simultaneously undergoing several phase transitions before 
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eventually transforming into HAP. 

This situation is seemingly incompatible with life as the persistent supersaturation of 
calcium and phosphate in biological fluids dictates unyielding sedimentation, at least in the 
absence of regulatory inhibition. Fortunately a regulatory mechanism does exist. The main 
machinery for preventing calcium phosphate sedimentation is the plasma protein fetuin- 
A (FA) [5], I2U1 - E3] . FA is an acidic protein that is found abundantly in blood as well as 
throughout all extracellular compartments. Maintenance of adequate levels of FA protein 
has been shown to be necessary for inhibition of calcification in relation to many disorders m 
[25]. 

FA interacts with the calcium phosphate mineralization process in several different ways. 
It can directly bind calcium, with each molecule able to weakly and reversibly bind approx¬ 
imately 12-15 ions [26]. Yet, this direct binding of calcium cannot be the main regulatory 
mechanism of FA as it would effectively reduce the supersaturation. There is also evidence 
that FA binds to pre-nucleation clusters [22], a possibility that has been analyzed [2Z], al¬ 
though somewhat contradictory evidence has also shown that the presence of FA does not 
affect the rate of nucleation of calcium phosphate clusters [2H]- Primarily, FA binds strongly 
to post-nuclear calcium-rich calcium-phosphate clusters, shielding them from further growth 
and imparting upon them enhanced colloidal stability so that they do not sediment. 

In this manuscript, we adapt mean-field classical nucleation theory (CNT) to look at the 
inhibition of mineral cluster growth by FA. We provide a quantitative description of the 
overall regulatory process and examine conditions necessary for stability. 


II. QUANTITATIVE METHODS 

The problem of understanding the combined process of mineralization and FA-induced 
inhibition is an example of a nucleation problem. Our approach to this problem is to 
use ideas from mean-held classical nucleation theory (CNT). In particular, we utilize a 
continuum approximation to the kinetic theory whereby we frame our problem using a 
series of serially-coupled partial differential equations (PDF). It is notable, however, that 
theoretical treatments of the inherently high-dimensional stochastic problem of nucleation 
and aggregation also exist [29H3T]. The mean-held theoretic CNT approach is ultimately 
motivated by the behavior of such stochastic treatments. 


3 


Our overarching goal in this section is to understand the kinetics of the concentration 
prohle for mineral clusters as a function of size and interactions with FA. To this end, in 


Section II A, we hrst solve for the nucleation rate (formation rate of critically-sized clusters), 
which provides a boundary condition for our PDF problem. Then, in Section II B[ we derive 
an effective growth rate (^^(s)) for the mineral portion of formed mineral and protein-mineral 
clusters as a function of their hxed mineral surface area s and number of bound FA proteins 
n. Then, in Section II C[ for a hxed cluster conhguration, we derive the shielding rate under 
the assumption that it is governed by diffusion-limited kinetics. Finally, in Section IIP we 


tie together the various components of our theory (nucleation, growth, shielding) into an 
overarching continuum model. 

For the reader’s convenience, we have compiled a list of the mathematical symbols that 
we use throughout this manuscript into Table |T} 


A. Nucleation 

CNT explains the emergence and evolution of colloidal phases in solutions through the 
development of a simple thermodynamical picture. The key element of CNT is the assump¬ 
tion that the emergence of a new phase carries an energetic cost due to the creation of an 
interfacial surface. 

Consider a mineral particle consisting of an integer number m mineral subunits, hereby 
termed an “m-cluster.” We say that this particle has volume V = mv, where v is the 
effective volume of each subunit “monomer.” Assuming that this particle is spherical, it has 
surface area s = (36h^7r)^/^m^/^, and radius r = In this manuscript, we use both 

m and s to parameterize the size of mineral clusters (see Fig. [^. We also will refer to the 
conversion between these two parameterizations as m{s) and s{m). 

For an m-cluster, CNT assigns as per the capillary approximation the free energy 

spherical 

AG = ^^ + mA^ = 7 /s +(1) 

where A/i is the molecular free energy per monomer in the cluster relative to in the solution 
(in units ksT), v is the volume per monomer, 7 > 0 is the interfacial surface energy per 
unit area (in units k^T per square meter), and / is a geometric factor than can be adjusted 
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Symbol Description 

s Surface area of mineral phase of cluster 

m Number of mineral monomers in a cluster, where 

monomer refers to Posner’s cluster Cag(P 04 )g. 
r Radius 

V Volume 

n Number of shielding proteins attached to surface of 

cluster 

7 Interfacial free energy per unit surface area in units 

k^T per squared meter 

/ Geometric correction for surface free energy for 

non-spherical pre-critical states 

poo Concentration of mineral monomers 

ps Saturation concentration 

A/j. — log(poo/ps )5 Chemical free energy per mineral 

monomer in units ksT 
AG Gibbs free energy in units k^T 

Cn{s^t) Concentration of mineral clusters 

of surface area s shielded by n FA monomers 
s' Shielded surface area 

Sn Amount of surface shielded (s') when n 

FA monomers are bound 
SSn Sn — Sn -1 for n > 1 

s*,m* Critical cluster size at nucleation 
Sp , Trip Critical cluster surface area and monomer number 
at sedimentation 

D Diffusivity of mineral monomers 

Dfa Diffusivity of FA monomers 

k- Dissociation rate per unit surface area for mineral 
0OO Concentration of FA monomers 

f, s, V Mineral monomer radius, surface area, volume 

OJ Dy/^ps/k- 

a SvDirpoo 

/3 

e Thickness of shielding layer (of FA protein) 

A \/4^_Dfa 0 oo 

TABLE I: List of mathematical symbols used in the manuscript for easy reference 

to account for non-spherical growth in the pre-nucleation stage as well as size-dependent 
variations in the surface free energy |32]. For constant supersaturation, Afi < 0 so that 
AG —)■ —oo as s —)■ oo, thereby thermodynamically favoring the existence of large clusters. 
Yet, as shown in Fig.|^a), the state of pure-monomers (s = 0) is also a local minimum of this 
free energy. The emergence of clusters is governed by kinetic rather than thermodynamic 
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s parameterization m parameterization 


FIG. 1: Cluster size parameterizations. Surface area (s) and monomer count (m) parameterizations 
are used interchangeably for expressing the size of mineral clusters. Clusters are assumed to be spherical 
and composed of an integer number m of monomers (Posner’s clusters). The cluster as a whole has surface 
area s. To denote the conversion between these two parameterizations we use functions m{s) : s ^ m and 

s(m) : m ^ s. 

considerations as an energy barrier of AGcnt corresponding to the free energy of a critical 
clnster with size s* mnst be overcome. This barrier is overcome when a clnster reaches a 
size m = m*, s = s*. The steady-state mean-held rate at which clusters reach this size is 
exponential in the magnitude of the energy gap and is known as the Zeldovich rate 

j* = Kexp (-AGcrit) (2) 

where k is a constant with units of concentration per time [33]. While the presence of 
pre-nucleation clusters in the calcium phosphate system violates the assumptions of CNT, 
Habraken et ah [H] showed that CNT is still applicable with the use of some minor modihca- 
tions that result in the reduction of the effective energy gap. Hence, we will assume for the 
purposes of this manuscript that critically sized clusters of size s* are forming spontaneously 
at some rate in the form of Eq. 

In the blood and extracellular compartments, fluid is under constant exchange. For this 
reason, we will also assume that the supersaturation is constant, and hence that j* and s* 
are hxed, and study the growth of clusters after their nucleation. 
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(a) Nucleation 



(c) Regulated growth 
n. 
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(b) Unregulated growth 
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FIG. 2: Nucleation, growth, and shielding, (a) Classical nucleation of an initial critically-sized 
spherical mineral cluster of surface area s* (« 3nm^). The activation energy AGcrit (units of ksT) 
corresponds to the energy of the cluster of critical size (of units surface area). A flux j* of nucleating 
particles is generated by the system under the condition of supersaturation, (b) Unregulated growth. 


Upon nucleation, particles grow uncontrollably at a size dependent rate vo{s) as defined in Eq. 22 (c) 


The attachment of proteins (shielding) alters the growth rate. The growth rate of the mineral 
Vn(s) depends on the size s of the mineral phase as well as the number of attached FA monomers n. The 
attachment of an additional protein to a cluster shielded by n proteins shields an additional surface area of 
size s„_|_i — Sn- Completely-shielded particles, where the surface area s is less than the shielding capacity 

Sn, do not grow. 


B. Growth of the mineral phase 

An m-cluster (of surface area s{m)) may find itself caked by a number of proteins, which 
effectively shield a surface area s' < s{m). Our immediate goal is to compute an effective 
growth rate for this particle assuming its shielding its fixed. We will assume that each 
successive protein shields a maximal surface area 6sn = Sn — Sn-i- In other words, if n 
proteins are attached, then a total surface area of size s' = Sn is shielded from further free 
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monomer adsorption. 

Due to surface reactions, this particle experiences an instantaneous net flux of monomers 
into its structure 

J = fc+(s, s')pr — X {s — s') (3) 

absorption dissociation 

where pr is the concentration of free monomers at the surface, is the dissociation rate 
per unit surface area, and /c_|_(s,s') is the absorption rate which is dependent on m as well 
as the free surface area. 

To begin, we will eliminate the unknown physical parameter function k^{s,s') by using 
equilibrium considerations to relate it to the other physical parameter k-. There exists a 
critical monomer concentration pm at which an m-cluster is at equilibrium with its surround¬ 
ings so that 

k+pra-k_{s - s') =Q. (4) 

At the equilibrium concentration, the free energy gap between clusters of size m and m -|- 1 
also disappears so that, 

5G = Ap + (36h^7r)^'^^7((m -|- = 0, (5) 


where Ap = \og{ps/pm) is the chemical potential, and ps is the free monomer concentration 
at saturation (where the solution is in equilibrium with in infinitely-large cluster). Eq. 
implies that 


Pm = p.exp [(36u27r)^/^7((m + 1)^/^ - m^/^)] 
= Ps exp 


/ Vm 

= Poo exp <; Ap 




1 _ 1 


f \ m 
, Ap 1/3 


— + 0{m-^) 

bm 


( 6 ) 


Substitution of pm from Eq. [^into Eq. [^yields the expression for 

Ps exp[-Apy/(s*/s)//J 
Eq. substituted into Eq. allows us to write the flux of monomers into the mineral cluster 


(7) 


J fc_(s — s') 


Pr - ps exp 


A 


p^exp 


( 8 ) 
















as a function of the physical dissociation constant fc_, the free monomer saturation concen¬ 
tration at saturation ps, and the free monomer concentration at the surface of the mineral 
Pr, and 


/5 = 




(9) 


The concentration pr is found through conservation of monomer mass by flux matching as 
in Fig. in the quasi-steady diffusive limit, where the FA protein forms a shielding layer 
around the mineral of thickness e. At a distance x from the center of the mineral, outside 
of the shielding layer {x > r + e), the free monomer flux obeys Fick’s law 

(r -h (5)(r -1- e) 


J = Airx^DdxP = 4:TrD- 


6-e 


'{Poo Pe) 


= 4:'KD{r + e){poo - Pe) + O 


r + e 


( 10 ) 


where the second equality is obtained by integration of the concentration from x = r + e io 
X = r -f (5, where 5 ^ r is the thickness of the diffusion layer. At the surface, the flux is 
given by Eq. Invoking free monomer conservation, by equating Eq. § with Eq. while 
also assuming that e is small relative to the characteristic diffusion length [pr ~ Ps), allows 
us to solve for the monomer concentration at the mineral surface. 


pr ~ Ps- 


A'KD{r + e)poo + fc-(s - s') 


A'KD{r -f £)p 5 exp (j + k_{s — s') 


exp 




( 11 ) 


In Eq. one sees that as D —)■ cx), the concentration at r goes to poo, as expected. Plugging 
the concentration from Eq. H into Eq. yields the growth rate 

A'KvD{r -h e){po 


V = vJ k_{s — s') 


pse 


( 12 ) 


AirD^r + £)pse^/^ + k_(s — s') 

While clusters of any size can be shielded, entirely-shielded clusters below the size s = Si 
are not of our concern because they are inert (recall that the hrst protein shields a maximal 
surface area of size si). Hence, we only wish to hnd the shielded growth rate when s > si. 
Making the assumption that e is small relative to r, for r = yjsjAir > yj si/dvr, yields the 
volume growth rate in terms of s, 

1 - s')(l - p^e^/^/poo) 


V 


with constants 


■y/lbTT + (s — s') 

oj = D\/~^ps/k_, 


(13) 
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and 


OL = SvDt^Poo- 


(14) 


Note that in the limit as s — s' -C this growth rate is snrface-limited, whereas 

in the limit as s—s' 3> growth is diffnsion limited. The parameters oj and (5 dehne 

the length-scale nnder which snrface-limited effects are signihcant. Typically, in nncleation 
problems, the snrface-limited regime is ignored as it is nsnally only important when the 
particle is small. For the pnrposes of onr system, we take a cne from the stndy of Trebonx 
et al. [13], which showed that aggregation of PCs is highly favorable energetically, and assnme 
that aggregation of post-nuclear PCs is diffusion limited. However, the surface-limited effects 
can become signihcant at larger particle sizes for clusters that are nearly-completely shielded. 
In this regime, we note that ~ 1, allowing us to use the simplihed growth rate 


1 — s') 

i/T^cnsP^ -K (s - s')' 


(15) 


The surface area growth rate is related to the volume growth rate through the chain rule. 


Wn(s) 



a(g - gn) 
a;sP2 (s - Sn)' 


(16) 


C. Shielding by the protein phase 

In this section, we compute the FA shielding rate for a hxed cluster conhguration. The 
shielding of the mineral phase by FA can be understood in a manner similar to the growth 
of the mineral phase. The overall adsorption rate of FA monomers onto the surface results 
from a balance between the diffusive supply and the surface reactions. Assuming hrst that 
mineral clusters have less mobility than FA, and denoting the diffusivity of FA by Ufa, one 
may use similar reasoning as in the previous section to hnd the overall attachment rate of 
FA. As before, we may express the flux into the surface as a balance between two competing 
reactions through a conservation law 

JpA ~ ^on0r('S S ) kogS (1”^) 

where kon is the binding rate of FA to the mineral per unit free surface area per unit 
concentration, k^g is the dissociation rate of FA, (poo is the far-held heat bath concentration 
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FIG. 3: Shielded diffusion-limited growth. FA protein forms a diffusion barrier of height e and 
surface area s' = Sn, where n is the number of associated proteins, around the mineral cluster. The 
absorption rate per unit surface area of the growth units depends on the local concentration p^. of 
growth units at the surface of the particle. Dissociation also occurs at rate fc_ per unit surface area. 
Neither absorption nor dissociation occur in the shielded region (red). For e sufficiently small, pr ~ Pe- 
The concentration pr is then determined through conservation of flux. Outside of the diffusion layer (of 
thickness S ^ r), the concentration of growth units approaches poo- 


of FA, and (pr is the concentration at the snrface. By Eq. 1^, we may also write the diffnsive 
flnx 

Jfa ~ F>fa\/^(0oo - 0r), (18) 


which, throngh conservation, matches the reaction flnx of Eq. 17 Eqnating Eq. and 


Eq. 1^ allows ns to solve for the snrface concentration of FA, 

^ _ IlFA\/47rS0oo + fcoffS' 

fcon(s - s') + UfaV dvrs 


(19) 
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Back-substituting Eq. into Eq [T^ yields the overall rate 

, -Dfa0 OO \/47rs fcoffs' 


JpA ~ ^on('S 5 )■ 


Goff'S • 


( 20 ) 


fcon(s - s') + 

At this stage we make a few simplifying assumptions - namely that the overall shielding 
process is diffusion-limited in the regime of most interest (where 0 < s' < s < Sp). First 
we remind the reader that our continuous formulation is an approximation of an underlying 
discrete system. For this reason, we note that the unshielded area s — s' cannot become 
inhnitesimal. This fact allows us to make the assumption that the binding reaction is always 
sufficiently fast such that fcon(s — s') DpA^dvrs, for s < Sp. 

Finally, we will assume that the rate of detachments is negligible (irreversible binding 
of FA to mineral clusters). Altogether, these assumptions allow us to write the simpler 
diffusion-limited shielding rate rule 


JpA ~ DpaV 47rs0 0 


( 21 ) 


Strictly speaking, the parameters 0oo and poo contained in these expressions are them¬ 
selves dynamical variables. Their evolution can be determined through mass conservation, 
as all changes are due to the balance between supply and consumption. We are interested 
however in the biologically-relevant situation where calcihcation is a local phenomenon cou¬ 
pled to global auto-regulatory processes that maintain supersaturation. For instance, fluid 
present in a knee joint is continually replenished through interstitial flow. That is to say, 
we set (poo and poo constant and examine the conditions for the regulation of sedimentation 
in this regime. 


D. Overall continuum model 

Classical work by Landau, Lifshitz [3l], dehned an advection problem to quantitatively 
describe the evolution of the cluster concentrations as clusters grow due to monomer ab¬ 
sorption. This work has been extended throughout the years [3S], and recently united with 
nucleation IMIET], which is introduced as an effective boundary condition. We further ex¬ 
tend this prior work by incorporating the effects of shielding. In this continuum approach, 
one may describe the evolution in size of the concentration prohle of clusters using an ad¬ 
vection equation, where the cluster growth rate provides an effective “velocity” or drift. 
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Overall, the dynamic concentration Cn{s,t) of clusters of mineral surface area s associated 
with n FA monomers is described for all non-negative integers n > 0 for s > by the 
partial differential equations indexed by n, 


dcn{s,t) 

dt 


Vn(s) 

/-^-V 

d I a{s - Sn) . 
as L Ll>S^/^ + (s — Sn) -I 

(Cn(s,t) - Cn-l(s,t)) 

"-V-' 

diffusion limited 


( 22 ) 


where Vn(s) is surface growth rate dependent on the number of bound FA monomers, s' = Sn 
for n > 1 and s' = 0 for n = 0, 

A = V^Dfa^oo, (23) 


and for notational convenience we set c_i(s,t) = 0. The right-hand-side describes diffusion- 
limited shielding of the mineral particles by FA protein, which is assumed to have high 
affinity for the mineral phase. The solution domain for Eq. 22 is shown in Fig. [^c). For our 
purposes, we will assume that the system starts at a reference time t = to ai the the initial 
state 

{ 0 s > s* 

(24) 

c“(sd s = s„ 

where c“(s*) is an equilibrium concentration set by the nucleation process. 

Critically-sized clusters (of size s*) are assumed to be created at the Zeldovich rate j* of 
Eq. 1^ This creation rate is balanced with consumption due to growth and shielding. As in 
Farjoun and Neu [SB], this growth flux is expressed in terms of the non-dimensional rate of 
number-growth, V/v of Eq. 
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Invoking their balance argument leads to the constraint 


j* = lim 


f^FA^oo \^d7r5 

■y/iTrZlpooS^''^ 


+ 


|co(s,t) 


(25) 


-)- s 

Hence, the combined effects of nucleation and shielding impose an effective Dirichlet bound¬ 
ary condition 


c“(s*) = Co(s*,f) 


= j* 


+ 


1/2 I 

USj + S* 


(26) 
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III. RESULTS 


In this section we construct a solution to the system of partial differential equations de- 


hned in Eq. and the boundary conditions dehned in Eqs. [M] - [26| We proceed hrst by 
non-dimensionalization of the problem formulation. Then, using the method of characteris¬ 
tics, we derive a sequential relationship between the solutions of the system. Approximating 
the solution of the system of equations, we analyze the steady-state behavior of the system 
in the limit where inhibition is sufficiently strong, from which we compute the overall rate of 
protein consumption and propose a criterion for effective inhibition of calcihcation. Finally, 
we parameterize our equations using experimentally measured values found in the literature. 


A. Nondimensionalization 

We seek a convenient non-dimensionalization of our serial system of PDFs describing the 
shielded growth problem. We begin by normalizing the surface area s, which ranges between 


the critical nucleation surface area 

s* and another critical surface area 

Sp which represents 

the surface area at sedimentation. 

Using these constants, we write the 

non-dimensionalized 

size variable 

s — s* 
s =- 

(27) 


Sp 

critical cluster size s*. 

« s* 

(28) 


Sp s^ 

shielded surface area. 




, 

(29) 


Sp s^ 

shielding levels s' = s„. 

« Sn 

(30) 


Sp S^ 

and surface-limiting parameter u 

OJ 

CO — 

(31) 


^/Sp s^, 

Rescaling time 

. _ a{t - to) 

(32) 


Sp — s* 

results in the series of non-dimensional advection equations of asymptotically unit speed. 
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5Cr,, ^ d 


(5 “1“ 


di ds _C()(s + + (s + s* — Sn) 

with non-dimensional shielding constant 

X{sp — 


= -X^/s + ■% (Cn(s, t) - Cn-l{s, t)) 


X = 


a 


where the concentrations have been scaled by the nncleation boundary condition 

A _ Cn(s(s),t(t)) 

Cn{S, 1) — 00 /' ^ ’ 

SO that the concentration of critical clusters is hxed 

co(0, i) = 1. 


(33) 


(34) 


(36) 


(36) 


B. Characteristics of the PDE system 


The non-dimensionalized partial differential equations of Eq. [^can be solved by invoking 
the method of characteristics sequentially for each PDE. The solutions to the PDEs contain 
the characteristic curves described by the equations 


ds 


{s + s^- s') 


(37) 


d£ a;(s-1-s*)P2 + (s + _ s') ■ 

The s — i characteristics, as shown in Fig. originate from points {io,s{io)). They follow 
the relationship 


i — tg — S — s(to) T 2 :UJ 


y/s -|- s* — \J s(to) -|- 


-1- us' log 


a/s + s* - x/F \/ H^o) + S* + 


ys + s* +y^s(to) + s^ - 


(38) 


Particularly, for unshielded clusters (where s' = 0), the last line of Eq. 38 
Along these curves, the concentration varies as 


IS zero. 


dcn 

df 


u 


s -1- s* s' 


2\/F+^ \Cj\/s -f s* s -f s* — s'] 
- A\/5 + 5* {Cn{i) - Cn-l{s,i)) . 


; (^) 


(39) 
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s 



FIG. 4: Characteristics for advection problem. Some sample s — t characteristics for the 
nondimensionalized PDE problem with Sn — s^, = 0.25. For (a) n = 0, and characteristics emerge from 
s = s*. For (b) n > 1, the characteristics emerge from the curve {(s,tc(s))}, where tc is given in Eq. 
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For the purpose of solving these equations, it is advantageous to invoke the change-of- 
variables 


u = \/s + s* 

(40) 

5 = '§n 

(41) 

ds = 2udu, 

(42) 


to reparameterize the curves as 


dt 

— = 2u 
du 


uu 


+ 1 


> 0 


Vn ^ "Y 


(43) 


from where it is evident that the relationship between i and u is bijective. Hence, we may 
use M as a proxy for f, hnding that the concentration prohles along these curves vary with 
u as 


dcn 

du 


ui 


M + Sr. 


— Sn CjU + — Sr 


+ 2\u 


2^^ + CjU — Sr 




+ 2\u 


2^"^ + CjU — Sr 
u‘^ Srt. 


c„_i(s(m),£(m)). 


(44) 


With the aid of an integrating factor, Eq. 44 can be written in the exact differential form 


+ CjU — Sn 

= 2XU^{U^ - Sni^"^ 


exp 


Xu^ \ Cj + -u 


exp 


Xu^ { u + -u 


Cn-l{u,i{u)). 


(45) 


The equation for the unshielded particle concentration cq corresponds to the homogeneous 
problem (sq = 0; c_i = 0). Respecting the boundary condition invoked by nucleation, as 
well as the initial Cauchy data, yields the solution 

UJ + \/s + s* \/^ 


co(s,t) = H (i-ic(s)) 


\/s + s* Cj + \/s 


exp 


X 


As* (Cj + 


exp 


A(s + s*) (cD + |\/s + s*) 


(46) 


where H is the unit step function and 

ic{s) = s + 2 Cj ( -x/s + s* - Vs 


(47) 
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is analogous to a “first-passage-time” for the formation of size-s clusters. 

To solve for the subsequent concentrations, we take advantage of the Cauchy data by 
initializing all characteristic curves along the curve {tc{s), s) given by Eq. 47, thereby setting 


= 0 at the left endpoint. Hence, each point (s, i) such that s > s' = Sn, and i > ids) lies 
uniquely on a single curve originating from 


\J s{io) + s* — 


^®±^exp 
- \Js+st-V¥ ^ 




+ 1 


^^^±^exp 
Vs+s*—VW ^ 


t—s—2ili (— \/I7) 


- 1 


(48) 


As growth of the mineral phase occurs more quickly in the unshielded clusters than in 
the shielded clusters, the s — t characteristics propagate quickest in the unshielded clusters. 
Hence, the hierarchy supp(c„) C supp(c„_i) C ■ ■ ■ C supp(co) holds. In fact, the supports 
of all functions Cn are equal, as necessitated by the coupling dehned by the right-hand-side 


of Eq. 33 The creation of size-s clusters of shielding n is driven more by the shielding of 
“n — 1 clusters” rather than the growth of “n clusters.” We use this fact, along with the 
presence of an exponential term within the exact differential of the right hand side of Eq. 45 
to formulate the ansatz 


Cn(s, i) = H [i- ids)) gnis, i) exp 


—A(s -|- s*) y Cj -\- — y/s -|- s, 


exp 


where 


As* I ca -|- — V s 


(49) 


gn{u) = 2A 


+ CjU — Sn 

(u^ — 


X 


- Sn)^‘^""gn-dQ)dq- 


W s{to)+s 


(50) 


According to Eq. 46 


9o{u) = 


u + u 


U CO -\- \/ s 


(51) 


For solving gn{sd)^ the lower bound for the integral in Eq. 50 is taken from Eq. 48 For 
solving the next equation 5 f„+i(s,t), all instances of s,f reparameterized by the variable u 


using Eq. 48 
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The iterated integrals of Eq. can be solved numerically through standard quadrature 
methods. Here we hnd some properties of the solutions to these equations before exploring 
their steady-state behavior, which is of the most interest to us. 


First, there is the question of whether these equations are well-posed. For u near vs 


one can invoke L’Hopital’s rule on Eq. 50 to hnd that 


lim gn{u) = 

XoJSn 1 

So, the solutions are bounded on the left. Now, we seek to find pointwise bounds for the 
solution away from the left boundary (for u > We note that the lower bound of the 


integral term in Eq. given by Eq. approaches ySn as t —)■ cx). Since the integrand is 
non-negative, the solution is bounded from above by the steady state solution 


gn{u) = 2A 


U + CjU — Sn 
{v? — 

nu 

X / q\q^ 


\\CjSt 


' ^/s{io)+s 


gn-i{q)dq 


< 2A 


U + CjU — Sr 


X [ q^{q^ - Sn)^^"^gn-i{q)dq 


^9niu) 


(52) 


By repeated applications of the Cauchy-Schwarz inequality, one sees that the integral 


term in Eq. 52 satishes the inequalities 

[ q\q^-Sn)^"^-gn-i{q)dq 


< 


< 


- SnY^^'^^dq 


1/2 


{u^ - Sn) 


4Ac2^577,-|-1 


- 1/4 

'u^-sY 


8 


I 111(vTijii) 

1/4 


A\CjSn ~\~ 1 

This computation gives us the pointwise bound on (/“, 


I 111■ 


9Y{u) < 


2A 


-|- OJU — Sr, 


{AXCjSn + 




X 


- sY 


1/4 


I 11 • 
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Using the fact that g^{u) = 0(1), it is easy to see by induction that gn{u) is bounded and 
smooth (behaving locally like a polynomial with order controlled by n) for u > y/s^, where 


it is of-note that g^ is bounded also in the vicinity of y/s^. Hence, by Eq. 49, each solution 
of Cn{s,t) is Schwartz-class, for all t>0. 


C. Steady-state behavior 


Our interest is in long-term behavior of the system. Observe that the solutions of Eq. ^ 
contain an exponential multiplicative factor that represents regulatory shielding. This shield¬ 
ing is strong provided that the term in the exponential is large, which is the case when the 
inhibition ratio obeys 

A(l + s,)3/2 ^ 1. (53) 

In physical units, this criterion can be expressed succinctly in terms of the concentration of 
FA protein needed, 

^ Poo ~ Ps (KAy 

(Poo 3> , (54) 

Dfa rup 

where rrtp is critical number of Posner clusters in a pure-mineral cluster at sedimentation. 
Note that if this condition were not to hold then a signihcant number of clusters of sedimen¬ 
tation size would form. Sedimentation would then occur until exhaustion of supersaturated 
species. In our subsequent analysis, we will assume that this condition holds. 

Since the overall solution is tapered by the exponential term which goes as or as the 
volume, we are most interested in the behavior of g^{u) in the vicinity of In this limit, 
we use the binomial theorem to approximate the integrals of the general form, for a, 6 G M, 

"" - Sn)% 

(a+l)/2 


"s/ Sn 






X 


XLbSn-\-b 


1 + — 


dx 


(n^ - Sn) 


Ac*)s7t,-|-6-1- 1 




+ 0 


_ Q ^AajSn+^+2 J 


(55) 


Eq. 55 allows us to evaluate g'^ to the leading order 

A(\/si -|- Cj) |-^ 


9T{u) 


T v s* Acusi 4- 1 


u 


Si 


-\- Cju\ . 


(56) 
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Through an inductive argument, one finds that for n > 1, 



( 57 ) 


■51 XuSj + 1 


We retain both — Sn and uiu in this expression because it is unclear which term is large. 


As n increases however, the Cju term will begin to dominate. This fact implies that surface- 
limited effects arise for large clusters, contrary to the situation for most other nucleation 
problems. 

D. Rate of FA consumption 

Since we know the rate of protein association as a function of n (the number of bound 
proteins) and s (the surface area of the mineral phase), we can compute the total rate of FA 
consumption. Denote Rjk the cumulative rate of mineral consumption in shielding particles 
of j FA monomers with size s G {sk,Sk+i)- We derive this rate first for the shielding of 
unshielded clusters of size at most Si, Rqq. Returning back to an integer parameterization 
of the size, it is clear that the total rate of consumption of FA for these clusters follows 


m(si) 



(58) 


The sum in Eq. 58 can be approximated by a left-Riemann integral so that 


m(si) 




m=m{st) 



(59) 


After transformation from m back to s, one finds that 



(60) 


Generalizing this result, it is easy to see that 



( 61 ) 
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Decomposition of solution domain 
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FIG. 5: Steady-state cluster concentration c“(s) and domain decomposition for computing 
the protein consumption rate. Green shading (darker is more concentrated) corresponds to increased 
steady-state concentration of mineral clusters of the particular size (given by horizontal axis) with the given 
number of attached FA proteins (given by the vertical axis). In the asymptotic case of strong shielding 


(Eq. 541, an exponential decay of concentration is seen according to size. From this solution, an overall 
consumption rate of FA protein can be computed by summing over the attachment rates , where Rnj 
refers to the rate of protein consumption by clusters of n bound proteins of size Sj to Sj+i. 


The total consumption rate of FA protein follows 


at nucleation 

w OO OO 

R ~ Aa/^c“(s*) + ^ ^ ^ ^ Rnj 

n=0 j=n 


( 62 ) 
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where the first term represents FA consumed in the instantaneous shielding of critical clus¬ 


ters. Substituting Eq. 49 yields 


R ^ (s*) < 1 -h exp 


f 2 

- 1 + Q x/'S 

a \ 6 


X 


E 


n=0 ■' 


S9n Hs)) 
■y/lhh^S^TT 


exp 


As 


ds 

a 

V 3 J\ 

i 


(63) 


using the change of 

variables 


where u{s) = \/s/{sp — s*). We approximate the integrals in Eq. 


63 


e = 


As 

a 





which we invert to hnd that 


s 

de 


3a 




V2A 


A 


- (w + a/s) ds 


a 

A 

a 


( 3a 


1/3 


cu 


UaF + 
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(64) 


For n > 1, we evaluate the integrals in Eq. 63 
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We approximate the integral in Eq. 65 by a sum of incomplete Gamma functions r(a, x), 

s{s — Sn + UJ^/s) 
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3 a V 3 
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( 66 ) 


where the remainder term is exponentially small by rationale of Watson’s Lemma. For 
n = 0, we have 


VTbhAAr 


exp 
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/ 2 ^A' 
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(67) 


( 68 ) 


In the regime where A/a —)■ cx), no clusters of size greater than s* form, and the rate of 
FA consumption matches the nucleation rate for AGP. Gonversely, if clusters of size greater 
than Si form, the FA consumption rate is strictly greater than the AGP nucleation rate. 
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E. Parameterization 


As we have mentioned, normal physiological calcium concentrations exceed supersatura¬ 
tion relative to the most thermodynamically stable phase of calcium phosphate. In fact, as 
we shall see, they also exceed supersaturation relative to AGP. 

Normal serum free ionic Ca^’*' concentration varies between 1.2mM and 1.3mM [38] . 
and normal total serum phosphate concentration varies between 1.12mM and 1.45mM. At 
pH = 7.4, one hnds using the Henderson-Hasselbalch equation that the concentration of 
free PO 4 ” lies between 3.7 and 4.9 nM. Various studies have explored the solubility of AGP 
relative to concentrations of its constituent ions (Ga^^, PO^^). By empirical formula, AGP 
has the negative-logio-solubility piP^ = 3pGa-|- 2 pP 04 ~ 26 at 310K and pH= 7.4 [39] . 
Using this calculation, one may compute the supersaturation relative to AGP, 

Hence, the supersaturation ratio is approximately ~ 1.3 in normal conditions. 

Estimates for the molecular weight of FA range from 51 — 67kDa, while the usual serum 
concentrations of FA range from 0.5 — l.Og/L. Thus, FA ranges in concentration between 
7/rM—19/iM in normal situations. We assume that rup = (P(IO^), where it notable that a 
single FA protein shields approximately 10^ Posner clusters [22]. For context, a cluster of 
size m = 100 corresponds to a diameter of approximately 4nm, assuming hexagonal close 
packing and using v = 0.3nm^ [22|. A diameter of 4nm is similar in extent to the size of FA, 
which has been measured to have a hydrodynamical radius of 4.3nm [301 • We also note here 
that the radius of a PG is approximately 0.4nm, so the ratio of the diffusivities between a 
PG monomer and FA molecule {D/Dpa), is approximately 10. 

To estimate the concentration of PG monomers (poo), we rely on indirect evidence as pre¬ 
cise quantihcation of these clusters does not appear to have been performed in the literature. 
We note that a study by Ghughtai et al. [TT] found that in physiological conditions approx¬ 
imately 6 % of solution Ga^"*" is present in calcium-phosphate complexes. A separate study 
has found that approximately 50% of calcium phosphate complexes have size consistent with 
Posner’s cluster m- An AGP nucleation study using 2.5mM free Ga^"'' and ImM K 2 HPO 4 
found spherical clusters of approximately 30 — 80nm in diameter after Ihr [ 6 ]. Assuming 
that the 80nm size corresponds to a cluster that nucleated soon after f = 0 , one hnds that 
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the concentration of PCs is at least 2nM in their preparation. Using this value, we estimate 
the an equilibrium constant for the formation of PCs 


k 


eq 


poo 

lCa^+|31POj-]2’ 


(70) 


hnding that pfceq ~ —16. This compution estimates nanomolar-range concentrations for PCs 
in the physiological range that we dehned above. 

With these rough estimates in mind, we may approximate an “inhibition-ratio” from 
Eq. as 


Dpoo 

The concentration of post-nuclear clusters is exponential in this ratio, suggesting that large 
values for I would inhibit calcihcation. We turn now to assays of calcihcation inhibition in 
order to validate this computation. 

Heiss et al. [13] assessed the inhibition of sedimentation in a highly supersaturated solu¬ 
tion of 20mM Ca^"*" and 6mM Na 2 HP 04 in a closed system. For this system, Eq. [TO] provides 
an estimate of poo = 2.5 x 10“^M. At 20pM, corresponding to / = (P(IO^), FA was shown 
to inhibit sedimentation fully over a time interval of days. At 1.5pM, corresponding to 
I = (P(IO^), FA was seen to initially inhibit sedimentation, but only for a period of two 
hours. In the latter case, we expect two things to be occurring. First, because the exper¬ 
iments are conducted in closed systems, exhaustion of FA is occurring over the long time 
span. Second, inhibition is exponentially weaker than it is in physiological settings implying 
a proportionally quicker rate of sedimentation. 


IV. DISCUSSION 


A. Protein consumption rate and implications 


The quantity R sets a minimum replenishment rate for new FA protein in order to 


maintain a steady concentration of FA, and hence colloidal stability. As seen in Eq. 26 


the parameter A is present in the denominator of c“(s*). As a result, to the leading order, 
R increases as A decreases. Failure to maintain this replenishment rate leads to decrease in 
0OO, the concentration of FA. A drop in (p^o further decreases A, thereby further exacerbating 
the situation (the less FA available, the more that is needed). Effectively, in the regime where 
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calcium and phosphate concentrations remain supersaturated, the number of FA molecules 
required to buffer each nucleating mineral particle increases as the concentration of FA 
decreases. 

Thus, even a small destabilization in FA replenishment can feed-forward to avalanche 
into catastrophic calcium phosphate sedimentation. This observation explains the experi¬ 
mental hnding that serum FA is often signihcantly depressed in systems exhibiting ectopic 
calcihcation, yet plentiful in the sedimented plaques [Hj. 


B. Assumptions, Limitations, and Extensions 

In our theoretical treatment of this topic, we have made some key simplifying assump¬ 
tions. Correspondingly, we have also limited the scope of our formulation and results. We 
reiterate that we are primarily interested in the earliest stages of the mineralization process, 
immediately after nascent nuclei have overcome the kinetic barrier and progression is gov¬ 
erned by thermodynamic considerations. For this reason, we do not consider later phases of 
calcium phosphates, as well as their nucleation through heterogeneous nucleation involving 
ACP precursors [19]. It is notable, however, that the transformations of calcium phosphate 
from ACP to HA have been a rich topic of research, and FA protein is known to interact 
with these phases as well, just as it interacts with ACP |40j . 

Biologically, we also assume that the specihc structure of FA is important in two ways. 
First, its hydrodynamical mobility determines the rate at which it is able to interact with 
mineral clusters. Second, its precise biochemical structure determines its propensity for 
strong interactions with post-nuclear clusters. We have assumed this mechanism to be the 
primary mechanism for FA-based inhibition for several reasons. 

Mainly, this mechanism has been suggested in the experimental literature [28]. Alterna¬ 
tively, FA has been reported to be able to bind calcium ions directly [26] , and hypothesized 
to bind pre-nucleation clusters (Posner’s clusters) directly [22]. Yet, binding free calcium 
ions directly would not inhibit nucleation without decreasing supersaturation. Binding of 
pre-nucleation clusters would decrease the nucleation rate for calcium phosphates, running 
counter to experimental evidence [28] . 

Aside from the experimental evidence, there are apparent advantages to the solubilizing of 
post-nuclear particles as both an effective and efficient strategy for controlling sedimentation 
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in super-saturated systems. While binding of pre-nucleation clusters would inhibit nucleation 
by increasing the energy barrier, it would require more inhibitory proteins to work because 
pre-nucleation clusters form at quicker rates than post-nucleation clusters. Experimental 
and theoretical evidence has conhrmed the presence and stability of Posner’s clusters in 
biological solutions. In our formulation, we assume that pre-nucleation clusters exist in 
quasi-steady equilibrium with free calcium and phosphate. We also assume that FA-free 
calcium and FA-pre-nucleation cluster interactions are weak and reversible so that we may 
ignore them. 

It is this specihcity of mechanism that makes FA a potent calcihcation inhibitor. Uniquely, 
experiments have shown that FA is the crucial protein for in-vivo mineralization inhibition. 
While other macromolecules such as albumin are known to interact with calcium phosphates, 
their physical attributes make them insufficient for this task. 

Albumin acts as a buffering agent for calcium in blood, helping to maintain Ca^"*" concen¬ 
tration in an analogous manner to maintenance of free H"*" ion concentration by pH buffers. 
For purposes of this study, the main effect of albumin is in setting the far-held equilibrium 
concentration of Ca^"'' and hence mineral monomers. We are ignoring interactions between 
the regulatory FA protein and other plasma proteins such as albumin. Although lacking 
in intrinsic capability, albumin has been shown to enhance the inhibitory properties of FA 
protein [IS] , however, their main effect is in later-stage stabilization of complexes containing 
multiple protein-mineral clusters |43j . 

In this manuscript we have also ignored other possible contributing factors to the overall 
mineralization process including interactions with other ions such as sodium, chloride, mag¬ 
nesium, zinc, or While these ions have been shown to inhuence mineralization, 

their importance in the early stages of nucleation is unclear. 

We also have not considered secondary interactions between mineral-FA hybrids, or the 
formation of calciprotein polymers. One of the goals of the present study has been to 
determine the content of the individual calciprotein monomers (the protein-mineral hybrid 
complexes we study in this manuscript). Observations by Wald et al. jlH] have shown that the 
stability and size of these secondary structures varies with the concentration of FA present 
in the system. A possible cause for this effect is the variations in the mineral to protein 
ratio in the clusters that we form in our model. The understanding of these calciprotein 
monomers gained from this study should prove useful in better-understanding the kinetics 


behind the formation of calciprotein polymers as well as subsequent phase transitions. 

Finally, while we have made an attempt to parameterize our model based on quantitative 
results from past literature, our numerical estimates remain rough guides and further exper¬ 
imentation is likely needed more-accurately parameterize our model. In particular, as the 
precise mechanism for calcium phosphate nucleation becomes better known, it will become 
easier to characterize the precise nucleation barrier and nucleation rate of calcium phosphate 
in order to better-understand the FA consumption rate. We have assumed for instance that 
the Posner cluster is the pre-nucleation cluster for this system, but recent literature has 
been mixed with respect to this hypothesis. We would like to emphasize, however, that our 
formulation is independent of the precise nature of these clusters as the bulk of our analyses 
focus on post-nucleation events. 

Looking more broadly at our work, the methodology that we have developed in this 
manuscript has potential in explaining a variety of solubility problems throughout biology. 
As an example, the system of stabilization, transport, and clearance of lipid molecules by 
HDL and LDL bears striking resemblance to the calcium-phosphate-FA system that we have 
analyzed in this manuscript. The formation of protein-non-protein complexes or colloids is 
a widespread feature of the homeostasis of solutions in biology. 


V. SUMMARY 

In this manuscript we have utilized classical nucleation theory to provide a quantitative 
description of the growth of calcium phosphate nanoparticles interacting with a shielding 
protein. In contrast with other theoretical work on similar systems, we have not neglected 
possible surface-limiting regimes of the process. Our quantitative description of the process 
provides an estimate of the critical concentration of shielding protein necessary for stable 
long-term inhibition of calcification, as well as an estimate of the total rate that the protein 
is consumed. 

Critically, we have found that disruptions of the ability to maintain the concentration 
of FA leads to increased overall consumption of the protein, and hence, exhaustion and 
sedimentation. 
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